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DESCRIPTION  OF  DETECTION  PROBABILITY  MODEL  AND  PROGRAM 

INTRODUCTION 

A program  has  been  developed  in  Code  3110D  to  calculate  the  probability  that 
at  least  N pulses  in  M consecutive  pings  will  exceed  a threshold  as  a function  of 
the  single  ping  signal  to  noise  ratio,  the  cumulative  false  alarm  rate  and  M. 

Before  describing  the  multiple  ping  program  a brief  description  of  the  program 
for  the  calculation  of  single  ping  detection  probability  will  be  given  since  it 
is  a subroutine  of  the  multi-ping  program.  Listings  of  both  ETRAN*  programs  are 
included  in  the  appendix, 

SINGLE  PING  MODEL: 

The  single  ping  model  is  based  upon  the  application  of  the  Neyman-Pearson 
criterion  to  signal  of  unknown  phase  added  to  white  Gaussian  noise  (ref.  2,3,^, 

5,7 »9)-  The  Neyman-Pearson  criterion  places  an  upper  limit  on  the  probability 
of  a false  alarm  and  seeks  to  maximize  the  probability  of  detection  within  that 
limit.  This  criterion  is  in  wide  use  in  both  radar  and  sonar  work  due  to  the 
difficulty  in  estimating  prior  probabilities  and  costs  which  are  req.uired  by  most 
other  criteria. 

The  optimum  receiver  (i.e.,  matched  filter)  resulting  from  application  of  the 
Neyman-Pearson  criterion  to  a signal  with  uniform  random  phase  and  additive  white 
Gaussian  noise  is  an  envelope  detector  followed  by  a filter  matched  to  the  known 
envelope,  as  is  well  known  (ref  2,3,1+, 5).  Rice  and  others  have  derived  the  ex- 
pressions for  the  probability  of  exceeding  a set  threshold  with  noise  alone  and 
with  signal  plus  noise  for  this  case  (ref.  3,5).  The  density  function  for  the 
signal  plus  noise  is  given  by  Rice  as 

lo(av) 

where:  v = the  random  variable  of  output  voltage  normalized  by  the 

rms  noise  density  from  the  matched  filter 

a = the  voltage  signal  to  noise  ratio  at  the  input  to  the  detector, 
(signal  amplitude/rms  noise  voltage) 

lo^’’')  = 'the  modified  Bessel  function  of  zero  order  and  argument  z. 

*ETRAN  is  a scientific  prograiriming  language  developed  by  NEL  Code  3110D  used  in 
Sonar  Performance  Prediction.  Reference  (ll)  describes  this  language. 


p(v)  = V exp 


2^  2 
V +a 


The  signal  processing  gain  may  be  included  in  the  expression  by  letting 


a = (2E/N  ) ' 
o 

vith  E the  signal  energy  and  the  noise  power  density.  This  is  an  upper 

limit  on  the  processing  gain. 

The  Neyman-Pearson  criterion  says  to  pick  an  acceptable  upper  limit  upon 

the  probability  of  false  alarm  and  maximize  the  probability  of  detection  under 

that  constraint.  Let  P be  the  probability  of  false  alarm  and  T*  be  the 

2 1/2 

normalized  voltage  threshold,  T*  = [R  /N^]  , at  the  output  of  the  detector. 

Then  the  probability  that  noise  alone  will  exceed  the  threshold,  R,  is  given 
by  Rice's  expression  (ref  5)  integrated  from  T*  to  infinity  with  a = 0: 


= / 


V exp 


(O)  d V. 


Making  a change  of  variables  Y = ^ and  noting  the  1^(0)  = 1 gives 

P^g^  = exp  - (T*)^/2 

Thus,  specifying  P^^  determines  the  threshold,  T*,  independent  of  the  signal 
level  (T*  is  a function  of  the  noise  power  denisty). 

The  probability  of  detection  (i.e.,  the  probability  that  signal  plus  noise 
will  exceed  the  threshold)  is  given  by: 


2^  2 
V +a 


j 


(av)  d V. 


In  general  the  integration  cannot  be  performed  analytically  and  curves  — ^ 
giving  Pp  vs  a for  various  values  of  P^^  have  been  generated  by  J.I.  Marcvmi  rr^ 
of  the  Rand  Corporation  (ref  9).  - i 
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The  single  ping  probability  program  calculates,  the  probability  of  detection 
as  a function  of  the  probability  of  false  alarm  and  the  signal  to  noise  ratio. 
For  the  cases  of  either  the  signal  to  noise  ratio  is  much  larger  or  much  smaller 
than  the  normalized  threshold,  T*,  the  asymptotic  expressions  from  Helstrom 
(ref  3)  are  used: 

for  a > T*  + 2.5 

Pjj  « 1 - (T*/a)^^^  exp 

^ (T*-a) 


for  a < T*  - 2.5 

Pp  = (T*/a)^^^  exp 

/iiT  (T*-a) 

For  the  case  where  the  signal  to  noise  ratio  is  close  to  the  threshold  T* 
the  modified  Bessel  function  is  approximated  by  the  first  three  terms  of  its 
power  series 

I (z)  = — (l+-i+  1—  + ...) 

° 8z  2!(8z)‘^  3!(8z)^ 


MULTIPLE  PING  DETECTION  PROBABILITY 


With  the  background  of  the  single  ping  probability  of  detection  calculations 
as  a function  of  false  alarm  probability  and  signal  to  noise  ratio,  we  shall  now 
examine  a different  detection  criterion  based  on  M pings.  The  following 
assumptions  are  made  for  this  first  case: 

a.  The  pings  are  independent  of  one  another  (i.e.,  the  probability  of  two 
pings  in  a row  exceeding  tne  threshold  is  the  product  of  the  probability  that 
each  of  the  pings  will  exceed  the  threshold.) 

b.  The  average  signal  to  noise  ratio  is  constant  over  the  M pings. 

c.  The  detection  criterion  is  that  N or  more  of  the  possible  M pings 
exceed  the  threshold  without  regard  to  the  order  in  which  N occurs. 

If  the  probability  of  exceeding  the  threshold  on  any  one  ping  is  p^  and 
of  failing  is  (l-p^)  then  the  probability  of  exactly  N pulses  out  of  M exceed- 
ing the  threshold  in  a particular  order  is 

o - N,,  ,M-N 

^1  Ps  (^-Ps^ 

Since  the  order  of  occurence  is  not  important  the  probability  of  exactly 
N out  of  M in  any  order  is  just 

„ /m\  ^ /m\  N ,,  xM-N 

^2=  N Pl=Nj  Ps  <^-Ps^ 


where 


{:)  ‘ 


is  the  binomial  coefficient 


/m\  _ M! 

\N/  ' N!(M-N)! 


The  probability  of  N or  more  threshold  crossings  out  of  M is  the  sum  of  the 
Pg's  since  the  events  are  mutually  exclusive. 


Prob  {N  or  more  /M) 


-I  (”) 


Using  this  equation,  the  probability  of  detection  using  a N or  more  out 
of  M criterion  may  be  calculated  in  terms  of  the  single  ping  detection  prob- 
ability. Similarly  the  probability  that  noise  alone  will  satisfy  the  criterion 
(i.e.,  a false  alarm)  may  be  calculated  in  terms  of  the  single  ping  false  alarm 
probability. 


In  applying  the  Neyman-Pearson  criterion,  we  wish  to  specify  some  "N  out  of  M" 
false  alarm  probability  along  with  a single  ping  signal  to  noise  ratio  and  calculate 
the  "N  out  of  M"  detection  probability.  Since  the  equation  for  the  output  false 
alarm  probability,  F^,  in  terms  of  the  single  ping  false  alarm  rate,  F^,  cannot 
be  easily  inverted  an  iterative  solution  is  obtained  by  estimating  F^  using  the 
relation  between  the  binomial  distribution  and  the  incomplete  Beta  function, 

I^(a,3).  An  approximation  to  the  inverse  of  the  incomplete  Beta  function  is  given 
in  ref.  1,  page  9^5. 


F 

o 


i (i) 


F^  = Ip  (N,  M-N+1) 


M 

i 
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F 

s 


N 

N+(M-K+1)  e^'' 


w = 


y^(h+X)^^^ 


- (— 

y2(M-N+l)-l  2N-1 j 


-t^/ 


2 dt 


1 

2(M-N+1)-1 
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The  resulting  estimate  of  is  used  in  the  binomial  sum  to  obtain  F^. 

The  error,  further  reduced  using  a correction  equation  of  the  form. 


F = F 
s si  — 


(c) 


This  correction  cycle  may  be  iterated  to  obtain  the  required  accuracy. 

Two  corrections  will  normally  give  something  less  than  \%  error  in  the  cumulative 
false  alarm  probability.  The  resulting  single  ping  false  alarm  probability,  F^, 
is  used  along  with  the  single  ping  signal -to-noise  ratio  to  obtain  the  single 
ping  probability  of  detection  as  discussed  in  the  first  section.  The  may  then 

be  used  in  the  cumulative  binomial  to  obtain  the  probability  that  the  "N  out  of  M" 
criterion  is  satisfied  when  the  signal  is  present. 


'■b/h‘1  (l)  Ps' 


The  result  of  such  a N or  more  out  of  M criterion  is  given  in  figures 

through  Each  figure  shows  the  probability  of  detection  vs  the  ratio  of  K/M 

with  each  curve  representing  a particular  value  of  M between  2 and  12.  Each 

figure  is  for  a constant  value  of  the  single  ping  voltage  signal-to-noise  ratio 

shown  as  the  value  of  A.  The  output  probability  of  false  alarm  in  all  cases  is 
_7 

approximately  10 


The  point  to  be  made  in  examining  this  data  is  that  for  any  fixed  value  of  N 
and  signal  to  noise,  increasing  M results  in  an  increase  in  the  probability  of 
detection  as  would  be  expected  since  the  available  signal  energy  has  been  increased, 
but  more  importantly  there  is  a "best"  value  of  N for  each  M.  An  expression  for 
estimating  the  best  N as  a function  of  M was  given  by  Berkowitz  in  ref.  2_^ 


N = .15  / M 


For  values  of  M between  2 atid  12  the  above  expression  results  in  a number 
within  plus  or  minus  one  of  the  observed  best  value  from  the  curves.  The  reduction 


6 


in  the  probability  of  detection  from  the  peak  value  due  to  using  a N from  the 
above  expression  is  small  due  to  the  slope  of  the  curves  changing  slowly  near 
the  peak  values. 

i 

Figure  ^ presents  some  of  the  same  data  in  a different  fonn.  The  probability 
of  detection  is  plotted  vs  M and  as  a function  of  the  signal  to  noise  ratio.  In 
this  case  N was  chosen  as  the  "best"  value  for  each  M (i.e.,  N=3  for  M=5).  In 
going  across  the  figure  along  same  constant  probability  of  detection  the  trade  off 
between  single  ping  signal  to  noise  ratio  and  M is  observed.  (i.e.,  for 
M=l,  N=1  requires  A=3.)  In  going  up  the  figure  along  a particular  value  of  M 
and  N,  the  increase  in  probability  of  detection  as  a function  of  signal  to  noise 
is  obtained.  Finally  along  the  lines  of  constant  signal  to  noise  ratio  the  increase 
in  probability  of  detection  vs  M is  obtained. 

This  analysis  was  based  on  the  assumptions  of  a constant  signal  of  unknown 
phase  with  additive  white  Gaussian  noise.  A less  restrictive  case  of  a signal  with 
uniformly  distributed  phase  and  Gaussian  distributed  amplitude  in  additive  white 
Gaussian  noise  may  be  solved  using  the  same  equations  as  above  if  the  signal  and 
noise  are  assumed  independent  and  the  signal  fluctuations  are  assumed  independent 
from  pulse  to  pulse.  Due  to  the  independence  between  signal  and  noise  and  the 
assumed  Gaussian  distribution  the  probability  distribution  of  the  sirni  is  also 
Gaussian  with  mean  equal  to  the  signal  mean  (noise  mean  assumed  to  be  zero)  and 
variances  equal  to  the  sum  of  the  signal  variances  and  the  noise  variances.  This 
is  the  same  form  as  the  input  for  the  case  of  constant  signal  except  the  variance 
is  increased  during  the  time  the  signal  is  present.  Thus  all  previous  results  are 
carried  over  if  an  equivalent  noise  is  defined  as  the  square  root  of  the  sum  of  the 
variances  and  used  any  time  the  signal  is  present  (not  when  determining  the  prob- 
ability of  false  alarm).  The  assumption  of  a constant  mean  signal  to  noise  ratio 
during  the  M pings  may  be  relaxed  in  practice  if  the  variation  is  small  and  the 
average  value  of  the  M single  ping  signal  to  noise  ratios  is  used  in  the  derived 
equations.  It  is  hoped  that  even  less  restrictive  cases  (i.e.,  signal  correlation 
between  pings)  will  be  programmed  in  the  near  future. 
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